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ABSTRACT 

We present a single epoch of high signal-to-noise ratio spectropolarimetry of 
the Type la supernova (SN la) 2004S taken nine days after maximum light. The 
flux spectrum is normal, but with the additional presence of high-velocity (HV) 
line features in both Ca II and Fe II. The object shows continuum polarization 
at the 0.4% level in the red, a value which appears to be typical of SNe la. The 
continuum data are consistent with a ~10% global asphericity in an axisymmetric 
geometry. Unlike previous observations of other SNe la with HV features, the HV 
features in SN 2004S show no strong polarimetric signature, though this may be 
due to the timing of our observations. Instead, the object shows line polarization 
features (P < 0.5%) that are rotated with respect to the axis of symmetry of 
the continuum. The line features are visible in Si II, Fe II, and Ca II, and 
appear to be narrowly confined in velocity space just above the photosphere. 
These polarization features are a result of compositional inhomogeneities in the 
ejecta. They may represent newly synthesized elements whose clumpy spatial 
distribution within the ejecta is distinct from that of the globally aspherical 
ejecta whole. 

Subject headings: supernovae: individual (SN 2004S) — polarization 



Introduction 



Carbon-oxygen white dwarfs that accrete enough mass to approach the Chandrasekhar 
limit can explode as Type la supernovae (SNe la). The first successful models of the explo- 
sion process invoked parameterized one-dimensional (I D) flame fronts propa gating outward 
from the center of the star as a subsonic deflagration (INomoto et al.lll984l ). However, the 
explosion process is fundamentally three-dimensional (3D) in nature as the propagation of a 
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deflagration front is marked by large-scale fluid motions and mixing due to Rayleigh- Taylor 
instabilities. Recent 3D pure deflagration models are not as successful at reproducing the 
observations because the simulated explosions are underluminous, underener getic, inefficient 



in burning the available fuel, and result in e jecta that are too well-mixed (IReinecke et al. 
2002[lGamezo et al.l[2003l : iKozma et aliboosh . 



Some modelers interpret these results as requiring the burning front to make a transition 



from a subsonic deflagration to a supersonic detonation as it beco mes turbulent (IKhokhlov 



199l[) . though the details of how thi s might happe n are u nclear ( Hillebrandt fc Niemeyei 
2000l ). Simulations, such as those of iGamezo et al.l (120051 ) . demonstrate that models with 



different prescriptions for the evolution of the burning front produce quite different ejecta 
structures. Observational probes of the geometry and structure of the ejecta are therefore 
necessary in light of the lack of an a priori theoretical understanding of the explosion. 

The deviations from spherical symmetry seen in the simulations discussed above are a 
result of the ignition process and the random physics of turbulent flame propagation. More 
systematic departure s from sphericity could be a signature of the progeni tor systems, which 
are largely unknown (IBranch et al.l 1 19951 : iLivid 1200 ll : iNomoto et al.l 120031 ) , or their environ- 
ments. Rotational flattening of the progenitor white dwarf, interaction with circumstellar 
material, or t he process of merging two white dwa rfs could lead to SNe la with large-scale 
asphericities JHowell et al.lboOll : Iwang et al.ll2003h . 



Spectropolarimetry provides a unique opportunity to probe the geometry of SNe la. 
Outward streaming radiation in a supernova atmosphere acquires a polarization when it 
is scattered by electrons. Distant supernovae are unresolved by the observer, so the net 
polarization of the received light represents an integration of emitted flux over the visible 
surface of the object. In the case of a spherical supernova, symmetry demands that the total 
polarization be zero. Non-zero ne t polarization is therefore a sign of some deviation from 
spherical symmetry in the system (jShapiro &: Sutherland! Il982l ). 



Global asphericities of the photosphere, such as e l lipsoi dal deformations, are manifested 
as net continuum polarization (jjefferyl Il99ll : iHoflichI Il99ll ). Past investigations of SNe la 
have found typical intrinsic polarizations of order a few tenths of a percent in the contin- 
uum. The object with the highes t measured intrinsi c continuum polarization percentage is 
still t he subluminous SN 1999by (IHowell et al.ll200ll ). which showed at least 0.7% polariza- 



tion. 



Howell et al.l (l200ll ) found that the spectropolarimetric data could be modeled by an 



oblate spheroid with an axial length ratio of 1.17, seen nearly equator-on. The lower con 
tinuum polarizations measured in the half-doz en other SNe la w i th good data are consistent 
with global aspherici ties at the ~10% level (jWang et al.l l2003l . l2006l : [Leonard et al.l 12005 
Chornock et"aDl2006h . 
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Line polarization allows other approaches to the study of asphericity. Polarization mod- 
ulations and angle rotations a cross a spectral featu re could be a sign of an aspherical dis- 
tribution of that ionic species. [Leonard et al.l (120051 ) found that the angle rotations seen in 
spectral features of some SNe la could be due to clumps or ionization asymmetries in the 
ejecta. Even if the ejecta lack such compositional inhomogeneities, line polarization can give 
insight int o the global aspher icity of the ejecta. Photons scattered in a line are generally de- 
polarized (IHoflich et al.lll996l ). but in a purely ellipsoidal geometry we should expect strong, 
isolat ed lines to be associated with "inverted" P-Cygni polarization line profiles (jjeffery 
19891 ). The effect of many overlapping lines is to produce a decrease in the p olarization at 
the same angle as the continuum, as was seen in SN 1999by (iHowell et al.ll200ll ). These signa- 
tures of asphericity are largely independent of the uncertainties in the removal of interstellar 
polarization (ISP) that plague estimates of the continuum polarization. 

One exciting recent development in SN la research is the association of "high- velocity" 
(HV) features in the flux spectra with strong spectropolarimetric features. These HV fea- 
tures are so-named because they are signatures of a separate line-formation region at higher 
velocities than the normal photos pheric component. S eparate HV material was first noticed 



in the Fe II lines of SN 1994D bvlHatano et al.l (I1999D . but it may be ubiquitous in SNe la. 
Every object in the sample of iMazzali et al. I J2005bh showed evidence of HV absorption in 



the Ca II near-infrared (NIR) triplet if the spectra were taken at sufficiently early epochs. 

The HV Ca II NIR triplet feature in the otherwise-normal SN 2001el was responsi- 
ble for a 0.7% spectropolarimetric f eature that was rotated with respect to the continuum 
( Wang et al.ll2003 : Kasen et al. 2003). This feature was modeled as formi ng in a spatially sep 



arate region that did not share an axis of symmetry with th e continuum (jKaseri et al. 



Subs equent observations of the SNe la 2004dt and 2002bf (ILeonard et al 



2005 



12003). 



Wang et al. 



20061 ) also showed strong polarization modulations across the HV lines. The composition 
and origin of this HV material remain mysterious. Suggesti ons in the literatu re have in- 
cluded clumps of material newly synthes ized in the explosion (jWang et al.ll2003l ). signatures 



of a gravitationa lly confined detonation (Kasen fc 



stellar material (jWang et al.ll2003l : iGerardy et al 



'lewa 



2004 



2005), and interaction with circum- 



Mazzali et al.ll2005al ). In all cases. 



viable models for the origin of the HV material must recognize its intrinsically 3D nature, 
as revealed by the spectropolarimetric line features and angle rotations. 

Inverting the observed polarization data for one object to construct a model is highly 
underconstrained. A more fruitful approach is to collect data for a sample of objects to see 
which, if any, model can bes t match the range of observed polarimetric behaviors. Though 
observations are still sparse, [Leonard et al.l (120051 ) have shown that SNe la exhibit diverse 
polarimetric characteristics that may correlate with other photometric and spectroscopic 
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properties. For example, the two SNe la with the highest observed intrinsic continuum po- 
larizations are both subluminou s, and the objects with HV fea tures have shown the largest 
line polarizations. In addition, IWang. Baade. &: PatatI (120071 ) have claimed to find a cor- 
relation between the strength of the spectropolarimetric modulation across the Si II A6355 
feature and the light-curve shape (a luminosity indicator) for a sample of SNe la. 

In this paper we add to the small sample of SNe la observed spectropolarimetrically 
and present a single epoch of high signal-to-noise-ratio (S/N) observations of SN 2004S, a 
SN la with HV features that is a virtual twin of SN 2001el. The data were taken 9 d after 
maximum light and the observed continuum polarization provides evidence for asphericity 
in the ejecta. Deviations from axisymmetry, seen as spectropolarimetric modulations across 
prominent spectral lines, are the main result and focus of the discussion below. 



2. Observations 

SN 2004S in MCG -05-16-021 was discovered by R. Martin on 2004 February 3.542 (UT 
dates are used throughout thi s pape r; Martin & Briggs 2004) and spectroscopically identified 
as a SN la bv ISuntzeff et al.l (120041). Optical pho t ometr y of SN 2004S has been published 



by both Misra et al. ( 20051 ) and Krisciunas et al. ( 2007 ). with the latter authors adding a 
significant near-IR light curve. They both find SN 2004S to be photo metrically normal 



and w e adopt a date for B-han d maximum light of February 4.37 ± 0.25 (IKrisciunas et al. 



20071 1. IKrisciunas et al.l (l2007h emphasize that SN 2004S and SN 2001el were found to be 
photometric twins. In addition, the spectra were very similar, with both objects showing a 
strong HV component to the Ca II NIR triplet lines at early times. 

We observed SN 2004S on 2004 February 13.34 (+9 d after maximum light) for a total of 
3600 s using the polarimeter unil^of the Low Resolution Imaging Spectrometer (LRIS; Oke 
et al. 1995) on the Keck I 10-m telescope. The D560 dichroic beamsplitter sent blue light to 
a 400/3400 grism and transmitted red light to the 400/8500 grating. The 1.5"-wide slit gave 
a spectral resolution of about 9 A over the full observed wavelength range of 3040-9250 A. 
The slit was aligned at 180°, which was near the parallactic angle for the first observation. 
Conditions were good and the seeing was ~0.9" during these observations. 

Standard two-dimensional (2D) image-processing tasks were performed using the pack- 



^See http://alELinoana.keck.hawaii.edu/inst/lris/polariineter/maiiual/pol_v3.ps for the online 
polarimeter manual by M. Cohen (2005). 



- 5 - 



ages m as were the optimal extraction of spectra and wavelength calibration. We 

used our own routines for flux calibration and atmospheric absorption correction of the ID 
spectra (IMatheson et al.l |2000| ) . The high airmass (~1.6) ar id flxed position angle of the 
slit combined to make the effects of atmospheric dispersion (IFilippenkd Il982l ) apparent in 
our data. As the first observation was taken close to the parallactic angle, we adjusted the 
continuum slope of the combined fiux spectrum to that of the first observation by dividing 
the total spectrum by a low-order spline fit to the ratio of the two spectra. 



The spectropolarimetric reductions followed th e standard procedure outlined by lMiller et al, 

The zero point of the waveplate was 



( 119881 ) and implemented by [Leonard et al.l (1200 ll ) 
determined by setting the observed positioi i angle of the polar ization standard star BD 
+59 389 to the catalogued value of 98.09° ( ISchmidt et al.lll992l ). As a check on the po- 
larization angle, we also observed a lower-quality polarization standard, HD 127769. The 
literature measu rements for the polarization vary somewhat, and the published angles are 
54.0° and 52.7° (jClemens fc Tapialll990l : iMathewson fc Fordlll970l ). We measured an angle 
of 52.5°, which allows us to conclude that our observation of BD -t-59 389 gave a reliable zero 



point. Two null polarization standards, HP 1090 55 and HD 57702 (iBerdyugin et al.lll995 



Clemens fc Tapial Il990l : IMathewson fc FordI ll97oh . were also observed and had negligible 
polarization (P < 0.07%). 

Dichroic absorption by interstellar dust along the line of sight polarizes light from the 
supernova. We selected two stars, CD —31 3659 and HD 48649, to act as probes of the 
Galactic contribution to the total ISP and observed them immediately after the supernova. 
Both are A-type stars more than 150 pc above the Galactic plane and within 0.6° of SN 
2004S; they can be regarded as intrinsically unpolarized probes that are distant enough 
to sample al most all of t he dust along that line of sight out of our Galaxy, following the 
suggestion of iTranI (1l995l ). 



The same dust that polarizes also reddens the light, which allows us to anticipate 
the magnitude of the eff ect of ISP. The Ga lactic extinction along the line of sight has 
E(B — V) = 0.101 mag (jSchlegel et al.l Il998l ). Empirically, dust in our Galaxy shows a 
polarization efficiency {Pmax / E{B — V)) of typically 3% mag~^, with a maximum value 
of 9% mag~^ (Serkowski et al. 1975), so we should expect the polarization to be ~0.3% 
if the line of sight is typical, up to a maximum of 0.9%. We integrated the fiux-weighted 
Stokes parameters over 5050-5950 A to approximate the V^-band polarization {qv,uv) and 



^IRAF is distributed by the National Optical Astronomy Observatories (NOAO), which is operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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measured (0.03, 0.49) and (0.11, 0.30), respectively, for the two probe stars. The two stars 
have somewhat different polarization, so we averaged their Stokes parameters to form our 
estimate of the Galactic ISP. Smoothed versions of q and u were then subtracted from the 
SN 2004S observations. Uncertainties in the determination of the Galactic ISP are absorbed 
into the determination of the host-galaxy contribution to ISP below. 



3. Results 



3.1. Total-Flux Spectrum 

The high S /N total- flux spectrum is presented in Fi gure [1] along with a spectrum of the 
prototypical normal Type la SN 1989B (jWells et al.lll994l ). We have removed a radial velocity 
for the host galaxy of SN 2004S of 2808 ± 5 km s~^ (determined from 21-cm observations by 
Theureau et al.l (120051 )) from all plots in this paper. The SN 1989B spectrum is from 8 d after 
5-band maximum light and is very similar to our SN 2004S spe ctrum, confirming that the 
age inferred from the ph otometry is ac curate. IWells et al.l (119941 ) measured the photometric 
decline-rate parameter JphillipsI ^9931 1 for SN 1989B to be Am^^(B) = 1.3 1 ± 0.07 mag. 



compared to the value of 1.14 ±0.01 mag reported by iKrisciunas et al.l (120071 ) for SN 2004S. 



We shall use the fiducial SN 1989B spectrum as a reference to better understand 
SN 2004S. The most obvious difference between the two spectra is that the absorptions 
due to Ca II are significantly deeper in SN 2004S and extend to much highe r velocities. 
SN 2 004S otherwise appears to be spectroscopically quite normal for a SN la (IFihppenko 



(120071), 



19971 ). A more complete spectroscopic sequence is presented in IKrisciunas et al. 
along with comparisons to a suite of SNe la. 

Figure [2] shows a closer view of the spectral features due to Ca II and Si II. The Si II 
A6355 line in particular is the hallmark of SNe la and is frequently used to measure the 
photospheric velocity (e.g.. Branch et al. 1988). As shown in the middle panel of Figure [2l 
the absorption minimum of the Si II line in SN 1989B is at ~10,000 km s~^, a typical value 
for a SN la. The Si II feature in SN 2004S is noticeably weaker an d the minimum is at 
only 8900 km s~^, compared to the figure of 9300 km s~^ reported by ISuntzeff et al.l (120041 ) 
near 5-band maximum light. This velocity is actually quite low for a normal SN la at thi s 
epoch and is comparable to the lowest- velocity objects in the sample of iBenetti et al.l (120051 ). 
Only peculiar subluminou s objects such as SNe 1991bg and 200 2cx have shown lower Si II 
velocities at a similar age (ILeibundgut et al.lll993l : iLi et al.ll2003l ). 



By contrast. Figure [2] shows that the Ca II spectral features in SN 2004S were formed 
at higher velocities than in SN 1989B. The Ca II H&K AA3934, 3968 absorption feature in 
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SN 1989B has a minimum at 3795 A (11,600 km s~^) and a rather flat-bottomed absorption 
with the blue wing rising to a flux maximum at 3675 A (21,200 km s~^), where the velocities 
reported here are relative to the (^/-weighted line wavelength of 3945 A. In SN 2004S, the 
absorption is much deeper and the flux minimum is much bluer, at 3710 A (18,400 km s~^), 
with absorption extending outward to 3600 A (27,400 km s~^). 

The profile of the Ca II NIR triplet shown in Figure [2] is more complicated than the 
other two lines, with at least three distinct flux minima present near 8080, 8230, and 8385 
A. Two effects contribute to the complexity of this line profile. The first is that the fea- 
ture is intrinsically formed by a triplet and the wavelength separation of the components is 
relatively large. The minima at 8385 and 8230 A can then be identified with Ca II A8662 
and AA8498, 8542, respectively, with a blueshift of 10,000-11,000 km s~^. Absorption min- 
ima at similar velocities can be be seen in the reference spectrum of SN 1989B. SN 2004S 
shows an additional flux minimum near 18,000 km s^^ with no corresponding feature in 
SN 1989B, which we identify as the counterpart of the high-velocity absorption seen in Ca II 
H&K. There are therefore two Ca II line-forming regions in SN 2004S, one at the photo- 
sphere and a second one separated from the photosphere by 8,000 km s~^. No photospheric 
component is visible in the H&K feature, w hich is dominated by th e HV component. The 



spectral sequence of SN 2004S presented by iKrisciunas et al.l (120071 ) shows the evolution of 



the Ca II NIR triplet absorption profile from being dominated by HV material at early times 
to photospheric velocities at late times. 

The exact velocity range of this HV absorption is difficult to determine directly from 
the spectrum because of confusion caused by line blending. If the minimum at 8080 A is 
identified with A8542, then the apparent velocity is 16,700 km s~^, but at that velocity the 
A8662 line would produce an absorption at 8190 A, which is close a local maximum in the 
line profile. It is likely that the HV material is at slightly higher velocity (such as 18,000 
km s~^), but with a velocity spread sufficient to blend the two components of the line profile 
together. A spectral synthesis calculation is necessary to make a more quantitative estimate 
of the velocity spread of the HV material. The maximum velocity of the absorbing material 
is easier to isolate, as the blue edge of the feature is at 7880 A (24,000 km s~^ if identified 
with A8542). 

As with all supernova spectral features, blends with other lines due to the large Doppler 
shifts may affect the strength and apparent wavelength of any given spectral feature. In 
particular. Si II A3858 and O I A8446 are candidates for making a potential contribution 
to the H&K and NIR triplet features, respectively. We do feel confident in identifying the 
presence of a high-velocity absorption component in Ca II, especially with the simultaneous 
detection of a feature at a similar velocity in both the H&K and NIR triplet features. We 
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note that no HV component is obvious in the profile of the Si II A6355 hne. There are 
notches in the blue wing, but these are the early signs of incipient iron lines that strengthen 
with time, not HV features (see the spectral sequence in Krisciunas et al. 2007). 



3.2. Spectropolarimetry 



The last step of preparation before we can begin analysis of the spectropolarimetry is 
to correct the data for ISP in the hos t galaxy. Normal SNe la have well-defined colors as a 
function of age, so iMisra et al.l (120051 ) were able to estimate a total E{B — V) to SN 2004S 
of 0.18 ± 0.05 mag by comparison with templates. The numbers given above imply that the 
host galaxy contributes E(B — V) 0.08 mag of reddening ^ so we expect a typical host ISP 
of ~0.24% (ISerkowski et al.lll975l ). iKrisciunas et al.l (120071 ) used the long baseline of their 
optical through NIR observations to set a tighter constraint on the extinction. They derive 
a host extinction of Ay = 0.071 ± 0.043 mag, which would imply a host ISP of 0.07% for 
typical dust, up to a maximum of 0.21%. Unlike the Galactic ISP, we have no method to 
directly measure the host contribution, so we will have to appeal to some interpretation of 
the SN polarimetry itself to make an ISP estimate. This ISP estimate will not be unique 
because different assumptions about the intrinsic SN polarization can lead to different ISP 
estimates. 

Some of the SN 2004S spectropolarimetry data are plotted in the Stokes q — u plane 
in Figure [31 Each circle represents a wavelength bin 50 observed A wide and the size of 
each is proportional to the wavelength of the bin it represents, so larger circles are from 
longer wavelengths. Anticipating the results below, we only plot data points from sections 
of the spectrum that we believe are free from strong line polarization features. The exact 
wavelength choices for these "continuum" data points are marked in the middle panel of 
Figure [H The wavele ngth dependence of ISP is sufficiently weak at optical wavelengths 
( ISerkowski et al.lll975l ) that the effect of subtracting a small ISP is similar to choosing a 
new origin for the data in the q — u plane. The distribution of plotted wavelength bins is 
elongated, with a clear spread of at least 0.4% in polarization between the blue end of the 
spectrum (smaller points) and the red end (larger points), which must be due to intrinsic 
SN polarization. The dotted circle in Figure [3] has a radius of 0.21%, corresponding to the 
maximum expected host ISP. 

The line plotted in Figure [3] is a fit to the data points. The simplest decomposition 
of the polarization is to choose an ISP vector that lies on the line and to assign the re- 
maining polarization to intrinsic SN polarization that varies in amplitude as a function of 
wavelength, with most of the variation occurring along the axis defined by the line. The 
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direction indicated by the l ine could then rep resent the symmetry axis of the SN continuum 
as projected onto the sky (IWang et al.ll200ll ). This argument does not give a unique value 
for the ISP, so we will appeal to SN la theory and spectral modeling ca lculations to a ssume 
that the intrinsic polarization at blue wavelengths must be small. As iHowell et al.l (120011 ) 
have argued, many overlapping lines of iron-peak elements are the dominant form of opacity 
in SNe la at wavelengths below 6000 A. The opacity from these lines should destroy any 



polari zation caused by electron scattering in an aspherical SN atmosphere. iHowell et al. 



( I2OOII ) then cho se an ISP vector to make the polarization zero at the blue end of their data. 



Other authors fiLeonard et al.l l2005l : IWang et al.l l2006l : IChornock et al.l |2006| ) have chosen 



ISP vectors to make the intrinsic polarization of certain blue emission features zero. 

For this work, we will select an ISP vector {qisPiUisp) that lies on the line in Figure [31 
The exact point along this line is arbitrary, but we chose the point (—0.040, —0.159), which 
we have marked with a square, in order to make the polarization small at blue wavelengths. A 
host ISP of this magnitude is consistent with the redde ning constraints mention ed above. We 
subtracted an ISP with the wavelength dependence of ISerkowski et al.l (119751 ) and assumed 
Rv = Av/E{B -V) = 3.1. 

The formal polarization [P = (g^ + u^)^] is positive-definite and is therefore biased high 
for low values of the polarization. We can avoid this bias by continuing to represent the data 
by two Stokes parameters, but we are free to choose a more convenient parameterization 
of rotated Stokes parameters (RSP; Trammell et al. 1993). In this case, we will rotate the 
ISP-subtracted data to align one Stokes parameter (called qRSp) with the obvious axis seen 
in the continuum data in Figure [31 and marked by the line. This coordinate syst em is similar 



to the decomposition of polarization into a dominant and orthogonal axis by IWang et al. 



( l2003l ). except that we have explicitly excluded wavelengths of strong line polarization from 
the fit. 

The rotated SN 2004S spectropolarimetry is plotted in Figure HI After the rotation, most 
of the continuum polarization is in the qusp Stokes parameter. The polarization rises from 
near zero at 4000 A to a maximum of 0.4% at 6800 A before declining at longer wavelengths. 
Several broad undulations in the continuum polarization are apparent as well as narrower 
features associated with spectral lines seen in the flux spectrum. Five of these narrower 
features are also visible as strong negative departures from the continuum in the ursp 
parameter and hence as rotations in 6. Spectral windows away from these narrower features 
are mostly fiat and near zero in ursp despite the broad undulations in qnsp, justifying our 
decomposition of the polarization into these Stokes parameters. We now consider the line 
features in more detail. 
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3.2.1. Silicon 

We will first examine the polarization of the Si II A6355 line because it is the most 
isolated and can be used as a reference to understand the other spectropolarimetric features. 
As noted above, the continuum polarization of SN 2004S is confined to Qrsp, while Figure H] 
shows that the Si II line has excursions away from the continuum polarization in both Qrsp 
and ursp- We can immediately conclude that the Si II line-formation region does not share 
an axis of symmetry with the continuum. 

This point is reinforced in the upper-left panel of Figure [5l where data points from 
both the continuum and the line are plotted in the q — u plane. The continuum points are 
the same as in Figure [3l but rotated into the new coordinate system and shown in gray. 
The elongation of the continuum points along a line in the q — u plane was argued above 
to be representative of scattering in a geometry with a single well-defined axis. The data 
points in black with error bars are 50 A bins between wavelengths of 5900 and 6400 A. 
These line polarization points clearly do not share the same distribution in the q — u plane 
as the continuum. They are arranged along a linear axis that is rotated with respect to, and 
intersects the major axis of, the spread of continuum points. 

Another view of the same data can be seen in the leftmost panels of Figure El where the 
Stokes parameters are plotted as a function of velocity in the line in 20 A bins. The Si II 
A6355 line is located in the middle of a broad depression in qnsp, which has local maxima 
near 5700 and 6800 A. There is a narrower, sharp feature in qnsp located just blueward of 
the minimum of the line in the fiux spectrum. This feature is at a similar velocity as the 
largest negative depression in ursp- The fiux minimum of the line is at 8900 km s~^ and 
the greatest deviation of the line polarization from the continuum is in the narrow range of 
approximately 9,000-11,000 km s~\ just above the photosphere. 

Interestingly, as one moves to redder wavelengths from the fiux minimum of the A6355 
line, ursp increases from its minimum, crosses zero near 6350 A (~0 km s^^), and rises to a 
peak in the red wing of the emission component of the line. This positive excursion in ursp 
returns to zero (the continuum value) near 6700 A, corresponding to a redshift of ~16,000 
km s~^, which may be a spectropolarimetric indication of Si II forming in the far side of the 
ejecta almost as far out as the expansion velocity derived above for the HV line-formation 
region. 
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3.2.2. Calcium 

The next ionic species we examine is Ca II. The analysis of the flux spectrum above 
presented evidence for two Ca II line-formation regions, at 10,000 and 18,000 km s~^, making 
the shape of the NIR triplet in the flux spectrum complicated. The polarization profile of 
the NIR triplet shown in Figure H] is also complex. Two separate polarization features are 
seen in both Qrsp and ursp near 8165 and 8345 A. 

The presence of two line-formation regions makes interpretation of the polarimetric 
features ambiguous, but the rightmost panels of Figure [6] can clarify the situation. The three 
lines labeled "HV" in the top right panel mark the expected location of features associated 
with the three components of the triplet, blueshifted by the 18,000 km s~^ velocity of the HV 
line-formation region seen in the total-fiux spectrum. The stronger of the two polarimetric 
features is at the same velocity as the expected A8662 line, but no feature is seen at the 
expected position of the other two lines. A much more natural identification of the lines 
can be seen in the lower-right panel. The lines labeled "PV" mark the expected location of 
features associated with the NIR triplet if they occur at the same velocity (10,000 km s^^) 
as the strongest polarization feature in the Si II line. The redder of the two polarization 
features is formed over a velocity range of 9,500—12,500 km s~^ if it is identified as Ca II 
A8662. The bluer of the two features is then a blend of A8498 and A8542 forming at similar 
velocities. In this interpretation, the two polarimetric features are formed at velocities just 
above the local minima in the total-fiux spectrum (~11,000 km s~^), similar to the behavior 
of the Si II A6355 line discussed in the previous section. 

The data points plotted in the q — u plane in the lower-left panel of Figure are 50 A 
bins between 8050 and 8300 A and correspond to the polarization feature we associated with 
AA8498, 8542. Data points from the A8662 feature, spanning 8300 to 8500 A, are plotted in 
the lower-right panel. In each panel the data points are connected in order of wavelength by a 
line. The location of these two features in the q—u plane reveals complex spectropolarimetric 
behavior. Not only do these Ca II AA8498, 8542 data points show a spread of 0.4% in each 
Stokes parameter, which is itself indicative of deviations from simple axisymmetry, but the 
line connecting the data points shows that the path in the q — u plane is ordered in the form 
of a loop. This is in contrast with the Si II data, which form a linear feature. The A8662 
feature shows a similar 2D spread with a suggestion of a loop, though with lower statistical 
significance. 
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3.2.3. Iron 

Next we will consider the two features visible in Figure H] as deviations from the con- 
tinuum in both ursp and 9 near 4800 A. The spectra of SNe la in this wavelength region 
are formed by a complicated blend of lines (mostly Fe II, S II, and Si II), so the identity 
of any given spectral feature is ambiguous. The wavelengths of the strongest rotations in 
these features are at about 4745 and 4832 A. We first assume that the spectropolarimetric 
features are formed in a similar manner to the Si II line analyzed above. The narrow core 
of the polarization feature in the Si II line was centered at an expansion velocity of ~10,000 
km s~^. De-blueshifting the unknown features by this velocity implies rest wavelengths of 
4906 and 4996 A for our mystery lines. The Ca II NIR triplet polarization feature was formed 
at a velocity of 11,000 km s~^, which would imply rest wavelengths of 4922 and 5013 A for 
these lines. 

One reasonable identification is with Fe II AA4924, 5018 if those lines are blueshifted 
by ~1 1,000 km s~^. There are two major complications for this interpretation. The first 
is that those Fe II lines form a triplet with A5169. All three lines share the same lower 
atomic level, so they should all be present if one of them is. With this identification, the 
stronger spectropolarimetric feature is associated with A4924, the weaker one with A5018, 
and A5169 lacks a counterpart. This ordering by strength of spectropolarimetric features is 
the reverse of the ordering by weighted oscillator strengths (gf). Naively, one might expect 
the strongest line to produce the strongest polarimetric feature, much as the strongest of the 
three lines produces the strongest feature in the flux spectrum. 

The second complication can be seen in the flux spectrum in Figure [TJ A series of 
flux minima at 4698, 4772, and 4885 A are labeled as Fe II. If this identification is correct, 
the three components of the triplet are blueshifted by 14,000, 15,000, and 17,000 km s~^, 
respectively. These velocities are much greater than the photospheric velocity measured from 
the Si II line, similar to the HV Ca II features discussed above, and thus raise the specter of 
multiple Fe II line-formin g regions, one at pho tospheric velocity and one at high velocity. HV 



Fe II was first invoked by iHatano et al.l (119991 ) to model the early-time spectra of SN 1994D 



and has since been used by several auth ors to explain spect ra of SNe la before and near 



maximum light (e.g., Stehle et al. 2005). [Branch et al.l (|2006[ ) found that their spectral fits 



to a number of SNe la were improved with the addition of HV Fe II. Suggestively, they also 
found that the object in their sample whose fit required the highest optical depth in HV Fe II 
was SN 200 lei, which also had strong HV Ca II absorption. Spectral fits that lack HV Fe II 
have difficulty producing alternative identifications for these absorption minima in SNe la. 



These difficulties with the Fe II interpretation lead to the natural question of whether 
any other ionic species could be responsible. The first alternative identifications we consider 
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are Si II and Ca II because these species are known to produce polarized spectral features 
at other wavelengths in this object. Ca II has no strong spectral features in this wavelength 
region, but the Si II A5051 blend is a plausible candidate for creating a spectropolarimetric 



feature. IWang et al.l ((20061) identified a feature from this line in SN 2004dt, but in SN 2004S 
this single line would have to explain two polarimetric features. These features would be at 
expansion velocities of 18,700 and 13,300 km s~^, significantly different from the polarimetric 
behavior of the strong Si II line (A6355), so we reject the Si II interpretation. 

The most plausible remaining candidate ion is S II, which has a large number of tran- 
sitions in the optical. Blends of these lines, particularly between 4500 and 5700 A, make a 
contribution to the spectra of SNe la. The strongest spectroscopic feature in the flux spec- 
trum attributed to S II in SNe la is the "W" -shaped absorption near 5300 A (see Figure [I]). 
No strong spectropolarimetric feature is visible in Figure H] near this wavelength. An iden- 
tification of the polarimetric features with S II would have to explain why the weaker lines 
leave a polarimetric signature and the stronger ones do not. 

We will adopt the identification of these lines with Fe II for the rest of the paper, noting 
the caveats above. Data points from these two lines (50 A bins between 4650 and 5000 A) 
are plotted in the upper-right panel of Figure [51 The distribution of points in the q — u plane 
is elongated in a linear manner along a direction roughly parallel to the linear feature seen in 
Si II, possibly indicating a common origin for the spectropolarimetric features. The middle 
panels of Figure [HI show the same data in velocity space. The expected locations of the 
three members of the triplet if they were present at the 15,000 km s~^ velocity inferred from 
the flux spectrum are marked in the qRSP panel, while the ursp panel shows the expected 
position of features at 10,000 km s~^ as measured from the Si II polarization feature. The 
polarimetric line features do not fall at the same position as the HV notches in the flux 
spectrum, but they do correspond to the positions of two of the three lines at photospheric 
velocities. 

One possible explanation for the lack of a polarimetric counterpart to A5169 is to rely 
on the observed narrowness of the sharp spectropolarimetric feature seen in the core of the 
Si II line and at the photospheric velocity in the Ca II NIR triplet. Perhaps whatever 
non-axisymmetric asphericity (e.g., clumping) is the cause of the spectropolarimetric line 
features and rotations in the other species is narrowly confined in the ejecta. Stronger 
lines are optically thick through more of the ejecta, possibly including the less dense spaces 
between clumps. If the A5169 line (with the largest /-value) is optically thick over more of 
the ejecta than the other Fe II lines, then the signature of the asphericity could be hidden. 
Another possibility is that blending with another line (Si II A5051 at photospheric velocities 
is a possible candidate) destroys the polarization signature. 
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3.2.4- N on- detections 

It is also worth considering which hne features do not show spectropolarimetric features, 
as well as those that do. The S/N of our observations degrades rapidly blueward of 4000 
A due to declining efficiency of the spectrograph and increasing atmospheric absorption. 
Despite these mitigating effects, we are confident that Figure S] shows no sign of strong 
spectropolarimetric features associated with the Ca II H&K lines of the sort seen in the NIR 
triplet lines of the same species. One potential explanation is that, if our choice of ISP is 
correct, the continuum at such short wavelengths has negligible polarization due to the high 
opacity in lines from iron-peak elements. The additional opacity of the strong H&K lines 
then results in little additional depolarization. 

Two of the more interesting ionic species that lack strong polarimetric signatures are 
O I and S II. Oxygen in the ejecta of SNe la could either be unprocessed material from the 
progenitor white dwarf or the result of incomplete carbon burning. In ei ther case the oxygen 



should not share t he same spatial distribution as the silicon. SN 2004dt (ILeonard et al.ll2005 



Wang et al.ll2006l ) showed the combination of s trong polarization in the HV Si II line features 



and a lack of a feature at the O I line, which IWang et al.l (120061 ) interpreted to be a result 



of turbulent burning processes leaving clumpy intrusions of silicon in an oxygen-dominated 
substrate. The O I A7773 line in SN 2004S is of normal strength (see Figured]), with a flux 
minimum near a radial velocity of 9,000 km s~^. No strong polarization feature is visible in 
Figure m in urspi and there is marginal evidence for an increase in Qrsp- Unfortunately, this 
line falls near the strongest atmospheric absorption at optical wavelengths (the A band) and 
so the S/N is somewhat degraded. 

Sulfur, however, should have a spatial distribution similar to that of silicon, as both 
elements are formed by similar burning processes. S II is responsible for the characteristic 
"W"-shaped feature near 5300 A in SNe la before maximum light, as well as more minor 
features just blueward of that. At this late epoch, the "W" -shaped feature has started to 
fade and the redder absorption is being dominated by the increased effect of iron on the 
spectrum at these wavelengths. We see no strong modulation in qRSP at the position of the 
S II features and ursp definitely lacks the strong feature seen in the Si II line. Differences 
in polarization behavior between the two species despite similar spatial distributions could 
be a result of different optical depths in the lines, or the underlying continuum polarization 
could be modified differently by the blends of lines from iron-peak elements. 
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Discussion 



The most basic question to be answered by spectropolarimetry is whether or not an 
object is spherically symmetric. The strong wavelength dependence of the polarimetric data 
shown above is convincing evidence of intrinsic polarization and hence of deviations from 
spherical symmetry in the ejecta of SN 2004S. 

With the set of Stokes parameters that we have chosen to present the data, Qhsp would 
be an estimator of the total polarization if the system had a well-defined symmetry axis. The 
need for a second variable to describe the data, here the Stokes parameter ursp, is indicative 
of deviations from axisymmetry. The second parameter is significantly non-zero only for 
wavelength regions near certain spectral features. Equivalently, the scatter of the continuum 
data points about the best-fit line in Figure [3] is similar to that expected from the errors. 
This means the SN ejecta can be approximated as having a global, axisymmetric asphericity 
with deviations caused by a non-axisymmetric distribution of certain ionic species. 

Polarization in a SN atmosphere is induced by light scattering off of electrons. Thomson 
scattering is a wavelength-independent process, but our estimator of the total polarization, 
qRSP, shows a number of broad undulations with wavelength and a general decrease in 
the polarization level toward the blue end of the spectrum. This is due to the nature of 
the opacity in SNe la, which is dominated by a large n umber of bound-bound transitions, 
primarily of iron-peak el ements (IPinto Sz Eastmanll2000n . These transitions act to depolarize 
the emergent radiation (IHofiich et al.lll996l ). so the gradient of polarization percentage from 
red to blue reflects the increasing n umber of such dep olarizing transitions in the blue and 
ultraviolet portion of the spectrum (IHowell et al.ll200ll ). Note that this interpretation does 
depend on our choice for the ISP, which we will revisit at the end of this section. 

The detailed form of the undulations seen in qRSP is therefore a direct probe of the 
line opacity in the continuum-formation region of SNe la, as wavelengths of greater line 
op acity should show more d epola rization. Insp ection of plots of the expansion opacity given 



m 



Pinto fc EastmanI (120001 ) and iKasenI (120061 ) shows a local maximum in the line opacity 



near 5300 A and a minimum near 4600 A, similar to the wavelengths of local minima and 
maxima, respectively, in the polarization data shown in Figure HI This correspondence may 
only be suggestive, but can be verified by modern 3D radiative transfer calculations of the 
type described by Kasen et al. (2006). 

The amplitude of continuum polarization gives information about the global degree o f 
asphericity of the SN atmosphere (jShapiro fc Sutherlandlll982l : |Jefferylll99ll : lHofiichlll99ll ). 
SN 1999by was the first SN la with a firm detection of intrinsic polarization, at least 0.7% 
( IHowell et al.ll200ll ). The best-fit model of iHowell et al.l (12001! ) explained the polarization 
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as being due to an oblate ellipsoid with the major axis 17% longer than the minor axis 
and viewed nearly equator-on. Despite the sometimes large line polarizations, subsequent 
observations of SNe la have yielded detections of intrinsic continu um polarization that are 
smaller than in SN 1 999by, with values of 0.3-0.4% being typical JWang et al.l boosl . I2OO6 



Leonard et al.l l2005l : IChornock et al.l l2006l ). Converting the continuum polarization to a 
geometry is model-dependent (e.g., Hoflich 1991), but th e lower polarizatio ns of most SNe la 
may be due to objects with asphericities closer to 10% (IWang et al.ll2003l ). SN 2004S is no 
exception, with a minimum of 0.4% intrinsic continuum polarization. 

The polarization angle rotations described above in line features of Si II, Ca II, and Fe II 
are evidence that the geometric distribution of these ionic species within the ejecta differs 
from the geometry of the continuum-enaittin g SN atmosphere. The pure ellipsoidal models 
of SN 1999by explored by lHowell et al.l (I2OOII ) had only radial variations in composition, but 
they were able to successfully model the spectropolarimetric feature associated with Si II 
present in their data. The additional opacity present in the strong line relative to nearby 
wavelengths resulted in additional depolarization at a constant polarization angle. Such a 
model for SN 2004S could possibly explain the dip in qnsp seen in our data centered on the 
Si II line, but not the feature in uptsp- In particular, the loop in the q — u plane for the 
feature associated with the Ca II NIR triplet cannot be formed in a simple axisymmetric 
model for the ejecta. 



The HV mat eri al in SN 20 1el a lso showed such & q — u loop (IWang et al.l 12003 
Kasen et al.l |2003|) . iKasen et al.l (120031 ) explored several different geometries for the HV 
ejecta and found that q — u loops were a generic consequence of systems whose components 
had different axes of symmetry. For example, they found that an ellipsoidal shell whose axis 
of symmetry was misaligned with respect to an ellipsoidal continuum region produced such a 
loop. They also found that models invoking clumps in the ejecta could produce q — u loops. 
Clumps can produce polarimetric signatures by selectively blocking and depolarizing part 
of the light from the underlying polarized continuum. These clumps might not be oriented 
along the symmetry axis of the system, and hence produce polarimetric angle rotations. 

In SN 2004S, however, the main line-polarization features are being formed at photo- 
spheric velocities. This can be seen most clearly in the Si II panels of Figure [HI While 
modulations of the polarization exist at all wavelengths, the wavelength bins with the great- 
est polarization deviations from that of the continuum are located just blueward of the 
P-Cygni absorption minimum. The main source of polarization could be clumps of Si II- 
absorbing material selectively blocking portions of the polarized photosphere and confined 
rather narrowly in velocity space above the photospheric velocity. Support for this picture 
comes from the two separate polarimetric features seen in the Ca II NIR triplet. Explana- 
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tions for the line polarization must prevent the polarization feature associated with A8662 
from blending with the polarization feature formed by AA8498, 8542, just 4000 km s~^ to 



the blue. Depolarizing clumps in the ejecta confined to velocities less than 
could be responsible. 



43,000 km s" 



Leonard et al.l (120051 ) described spectropolarimetric angle rotations in line features of 



Si II and Ca II in SN 2003du that are at least broadly similar to these results. They ascribed 
the rotations to either clumps in the ejecta or ionization asymmetries of the sort invoked by 
Chugail (I1992I ) to explain observations of SN 1987A. Clumps of radioactive ^^Ni formed in 
the explosion could produce local variations in ioni zation and e l ectro n density that would 
be manifested as polarization features. In addition, IWang et al.l (l2003l ) describe a potential 



"skin effect" as Ni clumps produce a tangential component to the radiation field whose 
effects would vary as a function of wavelength due to the wavelength dependence of the 
photospheric radius. However, any clumps of ^^Ni would only have a modest effect on the 
ionization state of SN la ejecta and the significant path length of gamma rays in the ejecta 
would tend to smooth out the effects of clumpy energy deposition from radioactive elements. 
Therefore, we do not favor the interpretation where polarimetric effects are directly due 
to the clumpy distribution of ^^Ni instead of the clumpy distribution of those ions directly 
responsible for the polarimetric features. 



Clumps of intermediate-mass elements could be detectable by other means. iThomas et al 



( 120021 ) examined the consequences of clumping on variation in line profiles. In their param- 
eterized models, the observed homogeneity of Si II line profiles in SNe la meant that the 
characteristic size scale of perturbations in composition must be of order 1% of the size of 
photosphere or smaller. Larger perturbations would result in variations of the line profile 
when observed from different orientations that are larger than the observed variation in nor- 
mal SNe la, though we note that in this particular object the Si II line is somewhat weaker 
than normal ( §3.1j Krisciunas et al. 2007 ). A single clump in front of the SN photosphere 
can result in large observed polarizations (IKasen et al.ll2003l ). so perhaps the small value of 
the polarization ra o dulat ions in the lines is a result of many small clumps, consistent with 
the iThomas et al.l (120021 ) results. Future numerical work will make possible quantitative 
calculations of the effects of clumps of different size scales on the polarization spectrum 
dKasen et allbooel ). 



Next we turn our attention to the HV lines. The evidence presented above shows that 
the line-polarization features in SN 2004S are associated with material near the photospheric 
velocity, not the HV line-forming region. A small polarization modulation is visible in ursp 
in the blue wing of the HV Ca II NIR triplet feature. A formal integration over 7875-8075 
A (blueshifts of 17,000-24,000 km s~^ relative to A8542) gives a non-zero measurement of 
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ursp = 0.103± 0.025%, which is significantly smaller than the photospheric feature. Naively, 
this is somewhat unexpected as SNe 2001el, 2002bf, and 2004dt all showed strong (0.7% 



< P < 2.5%) polarizat ion features in their HV features (IWang et al.l l2003l : [Leonard et al. 



20051 : IWang et al.ll2006l ). The answer may be that we caught SN 2004S at a special epoch, 
when the HV material had sufficient optical depth (or a large enough photospheric covering 
fraction) to imprint an absorption feature on the spectrum, but not enough to result in a 
significant polarization feature. 



The time series of spectropolarimetry presented in IWang et al.l ( 20031 ) showed that the 
HV polarization feature in SN 2001el was strong at early times and disappeared after max- 
imum light as the feature in the fiux spectrum also faded. Our data are most closely com- 
parable in time to their 2001 October 9 epoch (+10 d after maximum light), in which 
the HV fiux feature was still present, but fading, and the polar ization feature was ra ostly 
gone. SN 2003du provides another example of this phenomenon. iGerardy et al.l (120041 ) pre- 
sented spectra at ea rly times that showed a strong HV feature, but the data analyzed by 
Leonard et al.l (120051 ). taken 18 d after maximum light, show no HV absorption or associated 
polarization feature. 

A compar i son b etween the polarimetry of SNe 2001el and 2004S is very interesting. 
Leonard et al.l (120051 ) found that the spectropolarimetric properties of SNe la were diverse, 
but may be correlated with the photometric and spectroscopic properties. In this case, 
Krisciunas et al.l (120071 ) demonstrated that these two supernovae have almost identical light 
curves and spectra, so we would like to know whether the similarities extend to the polarime- 
try as well. Fortuitously, our single set of spectropolarimetry is at almost the sa me epoch 
(+9 d) as one of the observations of SN 2001el presented b y IWang et al. J2OO3I) (+10 d), 
though we do not have a corresponding time series of data. Wang et al. ( 20031 ) used a dif- 
ferent method to determine the ISP (see below) so our results are not directly comparable, 
but we can make qualitative comparisons to the data presented in their Figures 9, 10, and 
11. 

The continuum polarization of SN 2001el has a maximum amplitude of 0.3% and it 
decreases rather sharply by a few tenths of a percent blueward of ~5400 A, which is where 
SN 2004S also shows a change in continuum polarization. At earlier epochs, the Si II A6355 
line in SN 200 lei shows deviations from axisymmetry in the form of a g — m loop with 
an amplitude of 0.5%, similar to what we see in our SN 2004S data, though the feature 
in SN 2001el mostly faded by +10 d to the ~0.1% level. Some difference in the silicon 
distribution between the two SNe should b e expected from the sl ightly weaker and redder 
Si II line in the flux spectrum of SN 2004S (iKrisciunas et al.ll2007l ). 



As noted above, the HV Ca II NIR triplet polarization feature in SN 2001el largely 
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disappeared between maximum light and +10 d. The photospheric component of the NIR 
triplet in SN 2001el has a few tenths of a percent polarization at +10 d, although it lacks the 
two separate features seen in SN 2004S. SN 2001el also shows no obvious sign of the q — u loop 
at photospheric velocities apparent in the Ca II feature of SN 2004S. A polarization feature 
(0.2% in each Stok es parameter) in SN 2001el near 4800 A on day +10 is not mentioned by 



Wang et al.l (120031 ). but appears to correspond to the features in SN 2004S which we have 
associated with Fe II. In summary, SN 2001el has slightly less continuum polarization than 
SN 2004S and weaker line-polarization features at the epoch most comparable to our data, 
though the Si II line in SN 2001el at earlier epochs is similar to our data for SN 2004S. 

The interpretation of the results presented in this paper is affected by uncertainty in the 
determination of the ISP. Our choice of ISP has been guided by the belief that the intrinsic 



polari zation at wavelengths less than 5000 A should be minimal, as detailed by lHowell et al. 



(120011 ). Another plausible method for estimating the ISP is to obtain observations at a 



very late epoch, when the electron scatterin g optical depth sh ould be negligible and any 



measured polarization can be ascribed to ISP. IWang et al.l (120031 ) obtained such observations 



for SN 200 lei and determined a different value for the ISP than was measured from the first 



method (IKasen et al.l 120031 ). with no obvious way to reconcile the two possibilities. We lack 



late-time observations of SN 2004S, so we cannot perform this test. 

In this work we have chosen the first ISP method due to the simplicity of interpretation 
and the lack of a compelling argument for another ISP choice, but for the rest of this section 
we will examine the consequences of relaxing our assumptions. If we do not choose an 
ISP along the main axis shown by the continuum data points in Figure [3l then the angle of 
polarization of the continuum will vary with wavelength. One possible alternative ISP would 
be near the point {qRSP,UB,sp) of (0.0,-0.3) in Figure [5l chosen to minimize the polarization 
angle rotation through the Si II line as well as the magnitude of the polarization at the 
minimum of the P-Cygni absorption, but this ISP vector creates problems of interpretation 
as well. No choice for the ISP can eliminate the q — u loop shown by the Ca II NIR triplet. 
If we know polarization angle rotations must occur across a feature due to one intermediate- 
mass element, then perhaps we should expect them in others as well. In addition, the 
total host-galaxy ISP correction would be uncomfortably large given the low host-galaxy 
reddening. 

Potential ISP vectors other than our chosen one necessarily lead to non-zero polarization 
at wavelengths below 4000 A, which would be an important result. Implicit in our discussion 
of the continuum polarization is a construction whereby continuum radiation is polarized by 
electron scattering and then subsequently depolarized by bound-bound line-transport pro- 
cesses in the outer part of the SN atmosphere. This picture is backed by detailed numerical 
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simulations of radiatio n transport in SNe la. iPauldrach et al.l (Il996l ) analyzed the ID de- 
flagration model W7 (INomoto et al.l Il984j ) and found that at maximum light the optical 
depth to electron scattering was of order unity near a velocity of 5000 km s~^. The apparent 
"photospheric" velocities measured from the minima of P-Cygni line features are closer to 
10,000 km s~^ due the effects of pseudo-continuous opacity from many weak lines of iron-peak 
elements. 

At the later epoch of our SN 2004S observations (+9 d), the surface of optical depth 
unity to electron scattering should have receded even more deeply into the ejecta. In reality, 
the competition between polarizing electron scattering and depolarizing lines occurs over 
an extended range of r adii. The predominance of line opacity over Thomson scattering 
( iPinto fc Eastman! l2000l ) means that at wavelengths of high line opacity, the last scattering 
of photons occurs at larger radii in the ejecta (e.g., Hofiich et al. 1998), which limits the 
opportunity for those photons to be subsequently scattered by electrons and polarized as 
they exit the SN atmosphere. 

These arguments are based on ID properties of SN la simulations. Non-zero polarization 
at wavelengths of high line opacity could be a sign of clumps of ^^Ni mixed out into the 
ejecta. Such clumps could result in local enhancements of ionization and hence electrons 
that could produce polarization at larger radii than expected. Significant outward mixing 
of ^^Ni would be an interesting constraint on explosion models and would also have effects 
on the photometric rise time and overall light-curve shape. Choosing an ISP vector away 
from the line in Figure [3] also produces a slow rotatio n of the cont i nuum polarization angle 
with wavelength that is hard to understand, though IWang et al.l (120031 ) posit a potential 
"skin effect" . The time dependence of the polarization effects due to these hypothetical ^^Ni 
clumps is unclear, though obviously our single epoch of observations cannot provide any 
constraints by itself. 



5. Conclusions 

We have presented a single epoch of high S/N spectropolarimetry of the Type la 
SN 2004S, which had an intrinsic continuum polarization of 0.4%. This observation adds to 
a growing body of evidence that normal SNe la have global aspheri cities of order 10-20%. 
The data points from the continuum are consistent with an axisymmetric geometry such as 
an ellipsoid. Polarization angle modulations across spectral features due to Si II, Ca II, and 
Fe II are evidence that the distributions of some elements within the supernova ejecta differ 
from that of the continuum. One natural interpretation of these data is that we are seeing 
evidence for clumps of newly synthesized elements. These clumps could be a diagnostic of 
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turbulent burning during the explosion process (e.g., Gar nezo et al. 2005). Previous work 



has shown non-axisymmetric structures in HV material (iKasen et al.ll2003l ). but here the 



deviations from axisymmetry are present close to the photosphere. 

We argued that the clumps in the ejecta are sufficiently narrowly confined in velocity 
space to prevent the Ca II NIR triplet polarization features from being blended together. This 
conclusion could be tested with a time series of spectropolarimetric observations of a future 
SN la starting well before maximum light. At early times these clumps might be hidden 
because the lines are optically thick at higher velocities. As the absorption minima in the fiux 
spectra recede through the ejecta, the clumps could become visible polari metrically at later 



times. In particular, the spectroscopic sequence of SN 2004S presented by iKrisciunas et al. 



(120071 ) shows that at early times the flux minimum of the Ca II NIR triplet line was at 



higher velocity than the polarimetric features we see at this later epoch. Spectropolarimetric 



data of the similar SN 200 lei before maximum light (IWang et al.ll2003l ) show no features at 
photospheric velocities in the Ca II NIR triplet (which has a strong HV component), even 
though photospheric polarimetric features are visible in the Si II A6355 line (which does not). 

A time series of spectropolarimetric data would also be useful in understanding the HV 
material. SN 2004S is the first SN la observed polarimetrically with prominent HV features 
in the total-fiux spectrum that lacked counterparts in the polarization. We argued above 
that this was likely a consequence of the timing of our observations, as the HV polarization 
features in the sim ilar SN 2001el also disappeared a couple of weeks after maximum light 



( IWang et al.l l2003l ). The HV material must decrease in optical depth as it expands, so 
determining the epoch at which its effects disappear in both the total-flux spectrum and 
the polarization data has the potential to constrain the optical depth or covering fraction of 
the HV material, key parameters in any model. Alternatively, we could be viewing the HV 
material from a special orientation, such as along an axis of symmetry. 

The deviations from axisymmetry in SN 2004S are responsible for subtle variations in 
the polarization across the proflles of some spectral lines and would have been unobservable 
at lower S/N. Spectropolarimetric observations have the potential to be a powerful probe 
of SN la explosion models; however, quantitative conclusions will have to wait for detailed 
radiative transfer calculations of current 3D explosion models (e.g., Kasen et al. 2006). In 
addition, interpretation of the continuum observations is hampered by uncertainties in the 
removal of ISP, even for objects like SN 2004S suffering small extinction. Progress can be 
made with observations of a larger set of objects, preferably behind minimal dust columns. 
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Fig. 1.— To tal-flux spectra o f SN 2004S and the normal Type la SN 1989B from 1989 
February 15 (IWells et al.lll994l ). The numbers in parentheses are ages relative to 5-band 
maximum light. Some lines of polarimetric interest are labeled. The major difference between 
the two spectra is that SN 2004S clearly shows an extra high-velocity (HV) absorption 
component in both of the spectral features due to Ca II. See Figure [2] for details of these 
line features. The locations of the members of the Fe II triplet are marked at a blueshifted 
velocity of 16,000 km s~^. 
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Velocity (1000 km s"^) 



Fig. 2.— Prominent absorption lines in SN 2004S (solid) and SN 1989B (dashed). The 
velocity axis in each panel is computed relative to the 5f/-weighted line centroid for each 
blend (3945 A, 6355 A, and 8579 A, respectively). A separate multiplicative factor has been 
applied to the flux scale of the spectra in each panel to match the flux peaks of the two 
objects in order to facilitate comparison of the absorption features. In SN 1989B, all three 
hues have flux minima at similar velocities (10,000-12,000 km s~^), while in SN 2004S an 
additional absorption component is visible in both of the Ca II features at 18,000 km s~^, 
well above the photospheric velocity defined by the Si II line. 
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Fig. 3. — Continuum polarization data for SN 2004S in the q — u plane. The data have 
been corrected for the Galactic component of ISP, as determined from probe stars. Each 
circle represents a bin of width 50 observed A selected from wavelength regions free of strong 
line polarization. See Figure H] for the definition of these regions. The size of the circles 
is proportional to wavelength, so a gradient in polarization is evident between the shorter 
and longer wavelengths. The distribution of points is consistent with a single dominant axis 
of symmetry for the continuum. The line is a fit to the data points and the box marks a 
point along the line chosen to represent the host-galaxy contribution to the ISP. The dotted 
circle has a radius of the maximum expected ISP from host reddening arguments. After 
subtracting the chosen ISP vector, the data were rotated to align the new Stokes parameter 
qRSP with the axis of symmetry represented by the line. The cross at (0.35, —0.2) shows the 
median statistical error bars of the points in the plot. 
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Fig. 4. — Spectropolarimetry of SN 2004S. The polarization data plotted here have been 
corrected for ISP and binned to 50 observed A per bin for display purposes. In addition, the 
Stokes parameters have been rotated to a new coordinate system (g^sp, m_r5p), as described 
in the text. The total- flux spectrum is plotted in the background of the top panel to guide 
the eye. Rotations due to line features of Fe II, Si II, and Ca II are visible in ursp and 9. The 
horizontal lines overplotted in the middle panel represent sections of the spectrum that are 
free from strong line polarization features and are used to define "continuum" polarization 
data points in Figures [3] and [51 Note that the angle is ill-defined for wavelengths less than 
4200 A due to the polarization data points being close to the origin of this coordinate system. 
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Fig. 5. — Polarization of various spectral features in the q — u plane. The gray circles in the 
background of each panel represent the same wavelength bins from the continuum windows 
marked on Figure H] and shown in Figure [3l but rotated to the new coordinate system. The 
continuum points are clearly extended along a single line in the q — u plane, indicative of a 
single axis of symmetry. The points representing wavelength bins associated with spectral 
features due to each ion (shown in black with error bars) fall outside the spread of continuum 
points, implying that the lines do not have the same spatial distribution as the continuum. 
The Si II and Fe II features show similar polarization behavior that may reflect a common 
spatial origin. The Ca II data points are inconsistent with a simple single-axis symmetry 
because they show a spread in both qRSP ^^nd ursp- The feature associated with Ca II 
AA8498, 8542 appears to form a loop in the q — u plane. 
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Fig. 6. — Stokes parameters versus velocity for three species. The total-flux spectrum is 
overplotted in each panel to guide the eye. The hne velocities used for the coordinates are 
relative to the g'/'Weighted line wavelength for each ionic blend. Si II A6355 shows a broad 

polarization depression with a superposed narrow feature, present in both Qrsp ^-nd u^sPi 
and centered at 10,000 km s^^, just blueward of the photosphere. To aid in the identification 
of spectropolarimetric features, overplotted on the Fe II and Ca II panels are the expected 
wavelengths of each member of the respective triplets, blueshifted both by the photospheric- 
velocity (PV) and the high- velocity (HV) line-formation region for each ion. The polarimetric 
features are clearly associated with the PV for both ions, except for the lack of a feature 
associated with Fe II A5169. There are no strong polarimetric features associated with the 
HV line-formation region. See text for details. 



